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SUMMARY 


A mathematical model of the heating associated with leakage within the 
Space Shuttle Orbiter wing-el evon cove was used to analyze the aerothermal 
response of the structure. The analysis was applied to wind-tunnel results 
and to shuttle-entry conditions. Flow in the cove was modeled as developing 
flow between parallel plates. The mass flow entering the cove was assumed to 
be the portion of the external- flow boundary layer that is next to the wing 
surface ahead of the cove. A parametric study described the mechanism by which 
the energy of the ingested mass is distributed as a function of the Reynolds 
number, wall conductivity, surface emissivity, and cove height and length. In 
general, the rate at which energy was transferred into the cove interior was 
primarily determined by the combined effects of the convection, wall capaci- 
tance, and internal-radiation terms. 

Correlation with wind-tunnel results indicated that the cold-wall convective 
heating rates, wall temperature, and gas temperature of the cove were predicted 
by the analytical model for short exposure times. Predicted transient thermal 
response of the el evon cove subjected to shuttle-entry conditions indicated that 
the allowable leak area at the cove seal was 20 percent less than that previ- 
ously indicated by experiment. To predict more accurate structural tempera- 
tures, the techniques used herein to represent the convective heat transfer must 
be incoporated into a more detailed structural model which would more accurately 
represent the cove geometry and thermal mass distribution. 


INTRODUCTION 

Shuttle-type entry vehicles and hypersonic cruise aircraft experience dif- 
ferential pressure across the wings that can promote hot boundary- layer gas 
ingestion into the wing-elevon cove region as illustrated in figure 1 . The 
cross section of the shuttle wing shows gas leakage through the cove, past a 
rub seal, and into the interior cavity. Predicted external-surface tempera- 
tures at the shuttle wing-elevon junction can be as high as 1280 K (ref. 1). 
Obviously, the internal structure cannot withstand exposure to ingested gases 
at such elevated temperatures. Therefore, the gas ingestion must at least be 
minimized by cove seals such as the rub seal shown in figure 1. Since the 
design of an absolute seal is improbable, considerable effort has been made to 
determine the effects of cove leakage. 

A number of experimental investigators (refs. 2 to 5) have studied heating 
trends in the wing-elevon cove region. Experimental pressure and convective 
cold-wall heating-rate data, obtained by using small wedge-flap models (refs. 2 
and 3), indicate that pressure and heating decreased along the cove. Similar 
heating trends were obtained (ref. 4) from short-time tests (approximately 20 s) 
by using a full-scale model representing a section of the windward surface of 
the shuttle at the wing-elevon juncture. However, after an initial decrease 
at the entrance, the pressure remained constant along the cove up to the flow 


restriction representing the seal at the end of the cove. In another experimen- 
tal investigation (ref. 5), a full-scale model that was geometrically and struc- 
turally similar to the shuttle cove was subjected to a simulated shuttle-entry 
environment (up to 1 800 s) to determine the allowable seal gap. The correct 
thermal mass characteristics were considered important to simulate the hot-wall 
convective heating conditions and the internal radiation which was assumed to 
be a major contributor to the heat load of the cove interior. The temperature 
of the internal structure was monitored during each entry exposure at various 
seal-gap heights to determine the allowable gap size. 

Although the thermal response of a wing-elevon cove has been studied exper- 
imentally, a real need exists for adequate analytical tools that can be used 
to predict cove response under actual flight conditions, and some efforts have 
been made to meet this need. For example, in reference 6 a detailed flow anal- 
ysis of ingested gas is presented and the results of calculations are shown to 
be in reasonably good agreement with the experimental measurements of refer- 
ences 2 and 3 for steady-state pressures and heating rates. In this analysis, 
the fluid is a perfect gas with constant molecular weight and specific heat; 
and the one-dimensional flew analysis is combined with flat-plate heating cor- 
relations and nonsimilar boundary-layer computer solutions. The mean enthalpy 
of the ingested flow was obtained by an integration of the external boundary- 
layer flew. 

The present study was initiated to develop a convenient mathematical model 
that would be useful in predicting the thermal response of a wing-elevon cove 
during flight. The model described herein incorporates (1) a simplified one- 
dimensional flow analysis of the gas leaking through the complex cove geometry, 
and (2) the time-dependent thermal response of the structure. The analysis also 
accounts for convection from the ingested external boundary-layer gas, conduc- 
tion along the cove walls, and internal radiation of the cove walls. The valid- 
ity of the model is established herein by comparison with the experimental 
results of reference 4 and by the rational trends obtained from a parametric 
study. The model is then applied to predict the thermal response of the elevon 
cove to shuttle-entry conditions. The analysis presented herein is based in 
large part upon the analytical model developed in reference 7. 


SYMBOLS 

A c cove area, m 2 

A i cove leak area, m 2 

B cove wall thickness, m 

Cp specific heat at constant pressure, J/kg-K 

D cove gap height, m 

Dh cove hydraulic diameter, m 

d cove leak-slot height, m 
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F(z) radiation view factor between cove wall and environment at each end 

of cove 

h heat- transfer coefficient, W/m 2 -K 

K(z) radiation view-factor term between cove walls 

k thermal conductivity, W/m-K 

L cove length, m 

I cove leak-slot length, m 

m mass-flcw rate, kg/s 

N number of finite-difference segments 

N Nu Nusselt number 

N Pr Prandtl number 

p pressure. Pa 

Ap differential pressure. Pa 

Q heating rate, W 

q heating rate per unit area, W/m 2 

R free-stream Reynolds number, m -1 

P VD h 

R c cove Reynolds number, 

M 

T temperature, K 

t time, s 

V velocity, m/s 

x distance along cove length, m 

y partial boundary-layer thickness (fig. 13), m 

z dummy argument for view-factor equation 

a angle of attack, deg 

6 elevon deflection angle, deg 

e surface emissivity 
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y 


absolute viscx>sity, N-s/m 2 
p density, kg/m 2 

a Stef an-Boltzmann constant, W/m 2 -K 4 

4) distance along cove length of opposite wall, m 

Subscripts: 
c convection 

env radiation to cove wall from environment at cove ends 

ex cove-exit value 

g cove-leak gas 

i initial wall value 

in cove-entrance value 

k conduction 

m convective mass transport 

max maximum 

o environment at cove ends 

opp radiation to cove wall from opposite wall 

out radiation from wall 

RSI reusable surface insulation 

r net radiation 

rub rub tube 

s wing surface 

sw stored in wall 

w cove wall 

00 developed channel flow 
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ANALYSIS 


A simplified mathematical model has been developed for determining the 
aerothermal response of the Space Shuttle Orbiter cove formed at the wing-elevon 
junction. The ingested gas is assumed to be a perfect gas in equilibrium and 
the gas properties vary with temperature. The ingested gas flow is assumed to 
be a one-dimensional, hydraulically and thermally developing flow between two 
parallel straight plates; the effects of the curved-flow path on convection and 
radiation are assumed negligible. The flow properties are evaluated locally 
along the cove. The model is subdivided into a fluid model and a wall model 
as shown in figure 2. 


Fluid Model 

The fluid model, illustrated in figure 2 (a), includes the convective mass 
transport of the fluid and the convection between the fluid and the wall sur- 
faces, but the conduction through the fluid and the energy stored in the fluid 
are both neglected. Laminar flow is assumed with the cove of unit width. 

Convective mass transport .- The convective mass transport or net energy 
transferred along the cove in the flow direction is given by 


2m ~ 2m, in 2m, ex 


3T g (x,t) 


me 


Pr 9 


3x 


dx 


(D 


where 


m 


= PgVgD 


Surface convection .- The surface convection between the fluid and one wall 
is given by 


Q c = h (x, t) dx [T g (x, t) - T w (x, t) ] 


(2) 


where 


and 


Njj u (x, t) 


— N NU, oo "** 


RcNpr 

0.036 

(x + 6 *)/D h 


1 + 0.0011 


^Pr 

(x + 6')/D h 


(3) 


as obtained from reference 8. The parameter 5' is a number (0.01 for the pres- 
ent analysis) that assures a realistically large value of % u (x, t) at the cove 
entrance rather than an indeterminate value. The symbol % Uf0t> is the Nusselt 
number for fully developed flew between uniformly heated parallel plates and 
has a constant value of 8.23 (ref. 9). 

Energy balance .- A first-order differential equation is obtained frexn the 
summation of the energy transfer terms defined by equations (1) and (2) and is 
given by 


8T g (x, t) 

2h (X, t) [ Tg (X/ t ) — T w (X, t) ] + PgCp f gV gD =0 (4) 

9x 


The coefficient of 2 denotes that two walls exchange heat with the same ingested 
mass flow. The gas- temperature distribution Tg(x, t) is obtained from the 
nunerical solution of this equation for a given wall- temperature distribution 
T w (x,t) and for the boundary condition where the entrance gas temperature 
Tg ( 0/ t ) is specified or computed. 


Wall Model 

The wall model, illustrated in figure 2(b), includes the conduction along 
the wall, the net- radiation exchange at the wall surface, and the energy stored 
in the wall. The convection to the wall is also included and is the same as 
that used for the fluid model (eq. (2)). 

Conduction .- The net conduction across the segment dx, as discussed in 
reference 10, is given by 


Qk “ k W B 


9 2 T w (x,t) 


3x 2 


dx 


(5) 
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Radiation . - The radiation exchange between surfaces was analyzed by 
neglecting reflected radiation. The radiation leaving the wall surface dx 
is expressed as 


Qout = ea dx T w 4 (x, t) 


( 6 ) 


The radiation arriving at surface dx from the segments of the opposite wall 
is expressed as 


Qopp 


ecr dx f K (x— 4>) T w 4 (<l>,t) d$ 
J 0 


(7) 


where K (x— c}j ) d<j) is the view factor between segments d({> and dx and is 
defined by the following equation based on reference 11: 


K (z ) d(|> 


D 2 d<j> 

2 (D 2 + z 2 ) 3/2 


(8) 


where z is the axial distance between segments (z = x - <j>). The radiation 
arriving at dx from both ends of the cove (or channel) is given by the 
expression 


Qenv = ea dx [ p < x > T o 4 (0,t) + F(L-x) T D 4 (L,t)] 


(9) 


where the view factor is defined by the following equation based on reference 11: 


1 z 

2 2 (z 2 + D 2 ) 1//2 


( 10 ) 


The temperatures T o (0,t) and T 0 (L,t) are the physical-boundary temperatures 
at each end of the cove. For the present analysis, T 0 at the entrance was 
specified as the initial wall value. The T 0 at the exit was set equal to the 
exit wall temperature T w (L,t). The net-radiation exchange for a wall segment 
is given by the following sum: 


Qr - Qout “ Qopp Qenv 


(ID 
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Energy stored .- The energy stored in a wall segment is given by the follow- 
ing expression: 


Qsw 


Pw c p, w® 


(x, t) 

* it 


( 12 ) 


Energy balance .- Summing the appropriate energy transfer and storage terms 
and linearizing the radiation term by using the previous time step t-1 yield 
the following expression: 


3 2 T w (x,t) 

h(x, t) [Tg(x, t) - T w (x, t) ] - k w B - £0 

8x 2 


T w J (x, t-1) T w (x, t) 


- J* K(x-4>) T w 4 ( cf), t-1 ) d<£ - F (x) T o 4 (0,t-1) - F(L-x) T 0 4 (L f t-l) 

3T w (x, t) 


PyCp, w 3 


at 


( 13 ) 


where the initial wall- temperature distribution T w (x, 0) (designated T^) is 
specified and the boundary conditions at x = 0 and x = L are, respectively, 


3T W (0, t) 
3x 


= 0 


and 


3t w (L, t) 

= 0 

3x 


Values of T w (0, t) and T W (L, t) are calculated at each time step by using a 
finite-difference form of the governing equation and by using terms evaluated 
at the temperature of the previous time step. The value for T w (0, t) is cal- 
culated by using forward spatial differencing typically shown in reference 10. 
Thus, 
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h(0,t) [T g (0, t-1) - T w ( 0 , t— 1 ) ] 

k„B 

[T w (0, t-1 ) - T w { Ax, t-1 ) ] 

Ax 2 


- ea 


N 

T w 4 (0,t-1) - ^ K(0-4>j) T w 4 (4>j, t-T ) A<t>j 

j=l 


- F(0-0) T 4 (0, t-1) - F(O-L) T 4 (L, t-1) 


Pw c p, w B 
At 


[T w (0, t) - T w (0, t-1) ] 


( 14 ) 


whereas the value for T W (L, t) is calculated by using backward spatial differ- 
encing (ref. 10): 

h(L, t) [T g (L, t-1) - T W (L, t-1)] 
k„B 

[T W (L-Ax, t-1 ) - T w (L,t-l)] 

Ax 2 


- ea 


N 

T w 4 (L,t-l) - ^ K ( L — 4> j ) T w 4 ( (J> j , t— 1 ) AcJ) j 


j=l 


- F(L-O) T o 4 (0, t-1) - F(L-L) T 0 4 (L,t-l) 


p w c p, wB 

At 


[t w (Lf t) - T w (Lf t-1 ) ] 


( 15 ) 


An iterative process on each time step is used to solve for the gas- 
temperature distribution T g (x, t) and the wall- temperature distribution 
Tw (x,t). At each time step the T g (x,t) is obtained from equation (4) based 
on a given T w (x, t) (either assumed initially or obtained from the previous 
iteration step). This T g (x,t) is then used to solve for T w (x, t) by using 
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equation (13). The iteration process on each time step is repeated until the 
values for T g (x, t) from two successive iterations are within an allowable 
tolerance. 


PARAMETRIC STUDY 
Base-Line Conditions 

The mathematical model described in the previous section was used to demon- 
strate and characterize typical results and trends when basic parameters were 
varied. Selected base-line conditions, many of which are typical of the shuttle- 
elevon-rcove design, were used for this study. The cove wall material was the 
reusable surface insulation (RSI) and the gas-leakage medium was air. Other 
base-line values are listed as follows: 


Cove height, D, cm 

Cove length, L, cm 

Wall thickness, B, cm 

Wall emissivity, e 

Entrance gas temperature, T g (0, t), K . 
Initial wall temperature, T^ , K . . . 

Exposure time, t, s . 

Cove Reynolds number, Rq 

Leakage flow rate, m, g/s 


1.27 
38.1 
0.64 
0.8 
1 770 
295 
500 
50 
1.1 


Calculated Results for Base-Line Conditions 

Typical cove gas- and wall- temperature distributions calculated from the 
present mathematical model are presented in figure 3 for exposure times from 
0 to 500 s. For this example, the gas temperature at the entrance (x/L = 0 
in fig. 3(a)) was fixed at T g (0, t)/Ti = 6, and the initial wall- temperature 
distribution (t = 0 in fig. 3(b)) was specified at T w (x, 0)/T^ = 1. These 
curves describe a gradual penetration of energy into the cove with exposure 
time, which is a reasonable expectation. Thus, the initial gas- temperature 
distribution (t = 0 in fig. 3(a)) shows that the finite energy of the ingested 
gas is transferred to the cold cove walls within the first 30 percent of the 
cove length, but by t = 500 s the gas temperature at the cove exit has 

increased by a factor of 3. As indicated in figure 3(b), the wall temperature 

at the cove entrance reached radiation equilibrium first (t = 100 s). Then, 
the interior wall temperature approached equilibrium causing the location of 
the maximum temperature to move inward. The cooler wall temperature nearest the 

cove entrance can be explained as an effect of more direct radiation cooling to 

the external environment. However, convection appears to be the primary factor 
in control of the wall temperature through most of the cove, since the tempera- 
ture of the ingested gas was always higher than the wall temperature. The tem- 
perature distribution at t = 500 s was selected as the base-line condition for 
this parametric study to represent the total heat load of a shuttle entry. All 
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examples shown herein include all modes of heat transfer, but parameters which 
affect individual modes such as cove Reynolds number, thermal conductivity, 
emissivity, and cove gap height and length were varied to make significant com- 
parisons. The rate at which cove temperatures change is also a function of cove 
structural mass, but that parameter was not varied. 


Effect of Reynolds Numbers 

The effect of cove Reynolds number Rc on the calculated cove gas- and 
wall-temperature distributions at t = 500 s is shown in figure 4 for values 
of Rj, from 10 to 1000. As indicated, the cove gas and wall temperatures 
increased with R c . This result was expected since R c is primarily a measure 
of mass flow in the cove of constant cross-sectional area, and the amount of 
energy that enters the cove is directly proportional to mass flow when the 
entrance gas temperature is fixed. Furthermore, the total energy of the 
ingested mass represents the potential thermal load to the cove interior, and 
the ingested mass comes from a portion of the external boundary layer next to 
the wing surface at the cove entrance. Therefore, for constant external-flow 
conditions, an increase in mass flow into the cove (increased Rc) will also 
bring an increase in energy content since a greater portion of the external 
boundary layer containing higher energy flow will be introduced. The results 
of figure 4 also show that for each R c the wall-temperature distributions 
closely follow the gas- temperature distributions except near the cove entrance 
where direct radiation relief to the external environment caused lower wall 
temperatures. At Rc = 500 the calculations yielded gas and wall temperatures 
throughout most of the cove that were constant and close to the temperature of 
the gas at the cove entrance. Thus, an equilibrium condition was reached within 
500 s at those Reynolds numbers. At lower Reynolds numbers, cove temperatures 
would also be expected to reach equilibrium but at lower temperature levels and 
after longer exposure times. 


Effect of Wall Thermal Conductivity 

The effect of thermal conduction within the cove walls on the calculated 
cove gas- and wall- temperature distributions is demonstrated in figure 5 for 
values of k w /k WfR gj from 0.1 to 1000. The base-line distributions for shuttle 
RSI are given by the curves for •S//^w,RSI = ^ • When k w was reduced an order 
of magnitude to k w /k WfR gj = 0.1, the gas- and wall-temperature distributions 
were unaffected. However, when k w was increased, the results showed that the 
cove gas and wall temperatures toward the entrance decreased while these tempera- 
tures toward the exit increased. For an increase in k w by 3 orders- of magni- 
tude to k w /k WfR gj = 1000, the cove gas temperature became constant for 
x/L =0.3 and the cove wall temperature was constant over the entire cove 
length. These results were expected since large increases in k w allow faster 
penetration of thermal energy into the cove as energy absorbed by the walls is 
quickly distributed internally by increased conduction. For the shuttle the 
cove is constructed of RSI material; therefore, thermal conductivity of the 
walls should not adversely affect the cove thermal response. 
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Effect of Surface Emissivity 


The effect of surface emissivity on the calculated cove gas- and wall- 
temperature distributions is illustrated in figure 6 for values of £ from 0 
to 1.0. This range of values was chosen to cover the extremes of internal ther- 
mal radiation f rom zero at £ = 0 to maximum at £ = 1.0 where each cove sur- 
face absorbs and emits radiation at 100-percent efficiency. Reflected energy 
was neglected. These results indicate that by decreasing radiation (decreasing 
£) there was a corresponding reduction in radiation cooling that allowed the 
gas and wall temperatures toward the cove entrance to attain progressively 
higher values. (At £ = 0, the gas and wall temperatures at the entrance are 
equal. ) Also, the gas and wall temperatures toward the cove exit attained pro- 
gressively lower values as radiation was reduced. Thus, in the example, pene- 
tration of energy to the cove interior is retarded by reducing radiation. How- 
ever, as shown in figure 7, this favorable effect disappears as exposure time 
increases such that cove exit temperature increases for decreasing radiation. 

In figure 7, the cove-exit- temperature histories are plotted for various values 
of £ and over a time period that corresponds to a shuttle entry. These curves 
show that although reduced radiation is beneficial initially in retarding tem- 
perature rise in the oove interior, increased radiation beccmes advantageous in 
promoting radiation cooling during the latter part of the exposure when cove 
wall temperatures are high. The period of time over which radiation is benefi- 
cial would shorten for higher cove Reynolds number because less exposure time 
would be required to elevate the cove interior temperature. Therefore, the 
present analysis indicates that thermal radiation is an important factor in 
transferring energy from the ingested gas to the cove interior during a shuttle 
entry. A similar conclusion is indicated by seme earlier investigations. (See 
ref. 5.) 


Effect of Cove Height and Length 

The effect of varying the cove gap height D between the elevon and wing 
(fig. 2) on the calculated cove gas- and wall- temperature distributions is shown 
in figure 8. These results show that as cove gap height is reduced with the 
leakage rate constant, the cove gas and wall temperatures increase toward the 
entrance and decrease toward the exit in a fashion similar to that shown when 
the surface emissivity was varied. (See fig. 6.) As indicated by equation (3), 
reducing cove gap height yields higher heat- transfer coefficients. Conse- 
quently, a greater portion of the energy of the ingested gas is absorbed by the 
cove walls near the entrance as the cove height is reduced. The similarity of 
trends of these results to those of figure 6 suggest that transfer of energy by 
radiation to the cove interior also diminishes with decreasing cove gap height. 
However, as occurred with the effect of surface emissivity discussed in the pre- 
vious section, this favorable effect disappears and reverses as exposure time 
increases. 

The effect of varying the cove length L (fig. 2) on calculated gas- and 
wall- temperature responses with time at the cove exit as shown in figure 9 for 
cove lengths from 12.7 to 76.2 cm. Base-line length is 38.1 cm. These results 
show that the gas and wall temperatures at the cove exit decrease as cove length 
is increased. Although slightly lower heat-transfer coefficients are obtained 
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in the present analysis (eq. (3)) when cove length is increased, the primary 
factor responsible for the lower exit temperatures is the increase in thermal 
capacitance of the cove wall that accompanies an increase in cove length. Con- 
sequently, longer exposure times are required for penetration of ingested hot 
gas to the cove interior in coves with greater thermal capacitance. 


Summary of Parametric Study 

From the results of the parametric study it is seen that the model is able 
to account for the various modes of heat transfer and to delineate the mechanism 
by which the energy of the ingested mass is distributed as a function of all 
these modes. Thus, the present analysis shows that the rate at which energy 
is transferred to the cove interior primarily depends on convection, wall con- 
duction, internal radiation, and wall thermal capacitance. Although convection 
controls the cove wall temperature, energy penetration toward the cove exit can 
be retarded by low wall thermal conductivity, low internal radiation, and high 
thermal capacitance. However, when cove interior temperatures are high, radia- 
tion cools the cove interior. 


CORRELATION WITH TEST RESULTS 

The present mathematical model was used to analyze the results from aero- 
thermal tests of a large-scale model of a wing-elevon cove shown in figure 10 
and described in detail in reference 4. The test model was designed to obtain 
gas temperatures, pressures, and cold-wall heat-transfer rates within the cove, 
and, hence, exposure times were short. The general instrumentation distribu- 
tion is given in figure 11. Tests were made with a turbulent boundary layer 
on the wing pitched to an angle of attack of 1 2° and an elevon deflection angle 
of 10°. The free-stream Mach number was 6.9, the free-stream Reynolds number 
was 4.4 x 1 0 6 per meter, and the total temperature was 1888 K. Cove leakage 
was controlled by varying the size of a slot in the seal located at the exit 
of the cove. The leakage slot was either full or partial span as illustrated 
in figure 12. The table inset indicates the range of slot dimensions included 
in this analysis. 


Entrance Gas Temperature 

For proper correlation with the test results, the cove-entrance boundary 
conditions were modeled as described in figure 13. In figure 13(a) an illustra- 
tion of a simplified mathematical model of the present analysis represents the 
entrance boundary conditions including the ingested mass-flow rate m and the 
gas temperature Tg entering the cove. These entrance conditions are dependent 
upon the boundary-layer flow ahead of the cove and the external-wing surface 
temperature T s . The variation of T s as a function of time was computed by 
using wind-tunnel calibration data and the lumped-par ameter computer code 
described in reference 12. Next, the relationship between Tg and m was 
established for several values of T s by using a nonsimilar, integral method, 
boundary-layer analysis (ref. 13) of the external flow. The m within the cove 
was assumed to come from the portion of the external boundary layer of thickness 
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y (fig. 13(b)) next to the wing surface with no mixing at the cove entrance. 
The mass-related average temperature of the portion of the boundary layer next 
to the wing surface, which is greatly influenced by T s , was used as the 
ingested gas temperature at the cove entrance (fig. 13(c)). The entrance gas 
temperature, therefore, was a function of T s and m of the ingested flow as 
presented in figure 1 4. 

The solution used in the present analysis for correlating wind-tunnel test 
results incorporated iterative procedures that varied the entrance Tg and m 
until the relationship of figure 14 was satisfied at each time step. The solu- 
tion was started by assuming a value of the entrance Tg, and m was computed 
from the orifice-flow relationship by using measured cove exit pressure and an 
initially assumed exit temperature. The exit gas temperature and m were 
determined iteratively so that the exit temperature agreed with the cove gas 
temperature at the exit defined by the present analysis. For each assumed 
entrance T g , there were corresponding values of m and exit gas temperature 
that produced the measured pressure; and the corresponding value of T s was 
then used with figure 14 to obtain a new entrance T g to compare with that 
assumed. This iteration continued until the proper entrance Tg was estab- 
lished. Then this procedure just outlined was repeated for each time step to 
yield a transient solution of cove wall- and gas- temperature distributions. 


Cold-Wall Heating Rates 

Typical calculated cold-wall heating-rate distributions along the cove for 
various leak-slot areas are presented in figure 15 as solid-line curves and are 
compared with test data from reference 4. The ticked symbols represent the 
heating on the elevon side of the cove, and the unticked symbols represent heat- 
ing on the wing side of the cove. (See fig. 11.) The calculated heating-rate 
distributions agree with the data for the first half of the cove where the cove 
gap height is constant. However, the prediction is generally low for the cove 
transition region (varying cove height) and at the cove seal near the cove exit, 
since the analysis did not include the factors that would naturally produce 
increased heating associated with varying cross-sectional area. 

The cold-wall heating rates at three locations are correlated with mass 
flow in figure 16. The data include all of the leak slots tested that are 
listed in the table inset of figure 1 2. The correlation is best at the high 
flew rates where fewer experimental errors would be expected. 


Gas Temperature 

The experimentally determined cove gas temperatures obtained from exposed 
thermocouple beads in the investigation of reference 4 are presented in figure 17 
for various leak areas at two cove locations, x/L = 0.188 and 0.688. The solid- 
line curves are the predicted gas temperature obtained from the present analysis. 
The dashed-line curves represent the predicted thermocouple temperature calcu- 
lated by using a lumped-parameter model of the thermocouple and with the pre- 
dicted gas temperature. (The lumped- par ameter computer code is described in 
ref. 12.) The lumped-parameter model is illustrated in the inset of figure 17 
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and accounts for the mass of the thermocouple, the convection from the gas in 
the cove, and radiation losses frcm the thermocouple. For x/L = 0.188 
(fig. 17(a)), the thermocouple data are well below the predicted values from 
the present analysis. The predicted thermocouple response shows better agree- 
ment with the thermocouple data, which indicates that there was a considerable 
effect of thermal mass on the gas- temperature measurements. Moreover, the gas- 
temperature gradient along the cove near the entrance is very high so that the 
predicted gas temperature is very sensitive to the accuracy of x/L. At 
x/L = 0.688 (fig. 17(b)), the thermocouple data approached the predicted gas 
temperature (solid-line curves), and the predicted thermocouple temperature 
(dashed-line curves) accounted for the initial temperature response of the 
thermocouple. Therefore, the gas temperatures predicted by the present analysis 
seem reasonable. 


Wall Temperature 

The experimentally determined wall temperatures from reference 4 for vari- 
ous leakage areas at two cove locations, x/L = 0.06 and 0.78, are presented 
in figure 18. Wall temperatures were obtained from thermocouples attached to 
thin disks surrounded by thicker walls which produced local lateral conduction 
due to the local temperature differences of the wall. The predicted temperature 
for a wall thickness equal to the thin disk, shown by the solid-line curves, 
agreed with experiment for the first few seconds, but the experimental results 
then fell below the predictions because of lateral conduction. The predicted 
thermocouple temperature, shown by the dashed-line curves, was computed by using 
a lumped-parameter analysis (ref. 12) of these discrete locations (see inset 
sketch of wall thermocouple model in fig. 18(b)) to account for the conduction 
effect. The convective heat load was the same as that which produced the solid- 
line curves from the basic cove analysis. This adjustment to the present anal- 
ysis gives better agreement, but the basic analysis cannot accurately account 
for lateral conduction or variance in wall cross-sectional area. 

The elevon-cove thermal analysis appears to be useful in a general analysis 
of the cove heating problem. The prediction of wind-tunnel results from short 
exposure times indicates that the gas temperature entering the cove can be deter- 
mined from the external boundary-layer analysis. Also, the cold-wall heating 
rates and gas- temperature distributions of the cove can be reasonably predicted. 
The cove wall- temperature distributions are predictable if the wall structure 
can be represented as a plate with one-dimensional conduction. To predict more 
accurate structural temperatures, the techniques used herein to represent the 
convective heat load can be incorporated into a more detailed structural model 
which could more accurately represent the cove geometry and thermal-mass 
distribution. 


THERMAL RESPONSE OF GOVE TO SIMULATED TRAJECTORY 

The results frcm experiments on cove response to test conditions that simu- 
late a shuttle entry reported in reference 5 could be analyzed by the present 
analysis if the boundary-layer flow information for the test facility had been 
available. Unlike the procedure used for correlation of wind-tunnel test 


15 



results in the previous section, the procedure for predicting cove response to 
entry conditions includes a long-exposure condition where internal radiation 
would be a significant factor. Since the required boundary-layer flow informa- 
tion was not available, a comparison with reference 5 could not be made; con- 
sequently, the present analysis was instead applied to a shuttle-entry trajec- 
tory to demonstrate its computational capability to predict cove gas and wall 
temperatures. The results from this analysis are considered as a rough estimate 
of the cove thermal response because of the limited information available on 
local external-wing flew conditions and because of the inability of the present 
analysis to model the exact cove structure. Typical entry hot-wall heating-rate 
and pressure histories at the cove entrance from reference 5 are plotted in fig- 
ure 19. The wing-surface temperature was computed directly from the heating- 
rate history (assuming radiation equilibrium with an emissivity of 0.8), and 
in this case the wing temperature was used as the entrance gas temperature. 
Obviously, a more accurate gas-temperature value could tie obtained if knowledge 
of the approaching boundary layer were known and the methods outlined earlier 
in the analysis of the wind-tunnel results were used. The mass flow through the 
cove was computed by assuming a leakage gap area and the pressure drop across 
the oove seal for an orifice computation. The geometry, material properties, 
and gas properties were the same as the base-line values used in the parametric 
study. 

In figure 20, typical computed gas- and wall- temperature variations with 
time at the oove exit during entry are presented for various cove-seal leak-area 
ratios. The dashed-line curve in figure 20(a) represents the wing-surface tem- 
perature used as the cove-entrance gas temperature. From these curves, it was 
determined that as &i/A c was increased from 0.03 to 0.50, the maximum exit 
temperatures increased from 66 to 89 percent of the maximum entrance gas temper- 
ature. These results indicate that hot-gas leakage for even the smallest leak 
areas considered herein would produce a greater temperature at the cove exit 
than the current shuttle designs could withstand. However, proper account is 
not made for the actual thermal mass of the shuttle-cove interior structure. 

In considering the effect of structural mass using this analysis, addi- 
tional computations were made of the thermal effects on the shuttle cove-seal 
rub tube located at the oove exit. A lumped- parameter computation was made by 
assuming that the entire rub tube (illustrated by the inset in fig. 21) was at 
a uniform temperature resulting from a convective heat load over half of the 
surface area and from radiation away from the other half. The computed rub- 
tube temperature histories for various leakage areas and two initial tempera- 
tures are presented in figure 21 . The results in figure 21 (a) correspond to the 
exit temperatures of figure 20 for a "hot-start" where T^ = 295 K. The results 
in figure 21 (b) represent a "cold-start" condition where the shuttle is cooled 
to T^ = 1 85 K while in orbit. For each plot the dashed line represents the 
design rub-tube temperature limit of 450 K. Obviously, greater leak areas can 
be tolerated for the lower initial structural temperature (fig. 21(b)). The 
maximum rub-tube temperatures for these two initial conditions are plotted as 
a function of leak gap in figure 22. The allowable leak gap indicated for the 
"hot-start” condition is about 0.044 cm which is within about 20 percent of the 
experimentally determined value of reference 5. 
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The possible benefit of varying the cove internal radiation by means of 
cove gap height or surface emissivity, as suggested in the parametric study, 
was considered in the present application to entry conditions. For example, 
the cove exit wall- temperature histories for = 0.06 are presented for 

e = 0.8 and 0 in figure 23. Initially, the cove interior temperature was 
reduced by a reduction in radiation, but eventually this trend was reversed 
similar to that shown earlier by the parametric study in figure 7. In fig- 
ure 23 the maximum exit wall temperature, which occurred in the later portion 
of the entry, increased substantially as the internal radiation was reduced. 
The reduction of cove internal radiation would probably be beneficial for very 
small leakage areas (low mass-flow rates), lower than those included in this 
study. 


CONCLUDING REMARKS 

A mathematical model of the heating associated with leakage within a 
shuttle-type wing-elevon cove was used to analyze the aerothermal response of 
the structure. Flow in the cove was modeled as developing channel flow between 
parallel plates. The mass flow entering the cove was assumed to be the portion 
of the external-flow boundary layer that is next to the wing surface ahead of 
the cove. A parametric study described the mechanism by which the energy of 
the ingested mass is distributed as a function of the Reynolds number, wall 
conductivity, surface emissivity, and cove height and length. The results 
indicated that the rate at which energy was transferred into the cove interior 
was primarily determined by the relationship of the convection, wall capaci- 
tance, and internal radiation. Initially, low internal radiation reduced 
energy penetration into the cove; but when the interior temperature increased 
sufficiently, radiation helped cool the interior. 

Correlation of wind-tunnel results from another investigation indicated 
that the cold-wall heating rates, wall temperature, and gas temperature of the 
cove were predictable by using the present analytical model. Thus, assumptions 
used herein were correct. Predicted thermal response of the elevon cove sub- 
jected to shuttle-entry conditions indicated that the allowable leak area at 
the cove seal was 20 percent less than that previously indicated by experimen- 
tal results. To predict more accurate structural temperatures, the techniques 
used herein to represent the convective heat transfer must be incorporated into 
a more detailed structural model which would more accurately represent the cove 
geometry and thermal mass distribution. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23365 
August 11, 1 980 
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(a) Fluid model. 



(b) Wall model. 

Figure 2.- Illustration of math model of gas flow in cove. 
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Figure 11.- Cross section of model installed in 
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Figure 12.- Arrangement of cove seal and rub surface. Dimensions are given in centimeters. 
















Figure 16.- Correlation of cold-wall heating rates with mass flow. 
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(a) x/L = 0.06. 


T / T. 

W 1 



(b) x/L = 0.78. 

Figure 18.- Wall-temperature histories for various leakage areas 
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Figure 22.- Variation of maximum rub-tube temperature with leak gap. 
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Figure 23.- Cove-exit wall- temperature histories for various surface 
emissivities exposed to a shuttle-entry condition. Aj/Ac = 0.06. 
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